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Abstract 

Anisotropic spherically symmetric systems are studied in the connection and densi- 
tized triad variables used in loop quantum gravity. The material source is an anisotropic 
fluid, which is arguably the most commonly used source term in anisotropic stud- 
ies within general relativity. The gravitational-t-anisotropic fluid constraints are de- 
rived and analyzed and then quantum gravity inspired holonomy replacements are 
performed. The quantum properties of the fluid are dictated by the modified con- 
straint equations. Particular attention is paid to wormhole throats, as they provide 
a simplistic model of the structures thought to be ubiquitous in the quantum gravity 
space-time foam at high energy scales. In comparison to the purely classical theory, 
the quantum corrections act to increase the energy density of the fluid, which indicates 
that they may lessen the energy condition violation present in the classical theory. 
Related to this, in principle it would be possible to have scenarios where the classical 
solution yields everywhere negative (with a zero at the throat) fluid energy density but 
the corresponding quantum corrected theory possesses only small regions of negative 
energy density or even everywhere non-negative energy density. 
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I : INTRODUCTION 

In classical general relativity anisotropy has often played an important role from mod- 
els of stellar structure and cosmologies to more exotic solutions [l]-[6]. In the category 
of exotic solutions exhibiting spherical symmetry, anisotropy is almost always present in 
systems such as wormholes, gravastars, and even certain black hole models forming from 
gravitational collapse or singularity-free models. (See, for example, [7] - [H] and references 
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therein.) These sorts of anisotropic solutions usuahy involve strong gravity effects, and 
in some cases it may be argued that they are in the realm where effects from a quantum 
theory of gravitation should be considered. For example, the vicinity of the classical singu- 
larity formed during the gravitational collapse is almost certainly a domain where quantum 
gravity effects are expected to apply. More recently, Bojowald, Paily, Reyes and Tibrewala 
have performed a thorough, in-depth study of midisuperspace models containing horizons 
with scalar fields ^15j in order to study quantum effects in these space-times. 

In the arena of non-trivial topologies, a common rationale for studying wormholes is 
that they may provide a simple model for the space-time foam, a potential model for the 
vacuum in theories of quantum gravity, where not only the geometry but also the topology 
of space-time may fluctuate due to quantum effects, [16], [17], as illustrated in figure [H 
It has been argued that systems of wormholes may provide good models for such a foam 
|18j . |19j . [20j . Therefore, in the vein of wormhole-like structures, there is good reason to 
study these systems utilizing the variables of loop quantum gravity due to their relevance in 
topology change which may exist in the quantum gravity realm. In the context of black hole 
physics, topologically non-trivial horizons have already been studied within loop quantum 
gravity [21]- [21]. 




Figure 1: A possible picture for the space-time foam. Space-time that seems smooth on 
large scales (left) may actually be endowed with fluctuating topologies (repre- 
sented by handles on the right) due to quantum gravity effects. One of the 
simplest models for such a handle is the wormhole. 



There have been a number works produced studying the possibility of a change in the 
topology of space-time at the classical and semi-classical level [25j - |44j. It is known 
that if one allows for the possibility of topology change in space-time, there are arguments 
that if Vo and Vi are distinct compact 3-manifolds, there will exist a space-time whose 
boundary is comprised of the disjoint union of Vq and Vi [25] (see figure [2] for reference.) 
Regarding this, Geroch's theorem states that if a segment of space-time has Vq and Vi 
as boundaries possessing differing topology, then a singularity or closed time-like curves 
must exist somewhere on the manifold |26j . It has been convincingly argued that these 
singularities that arise in certain topology changing space-times are extremely mild [27J 
in the sense that the tetrad becomes degenerate (introducing metric degeneracy) but the 
resulting curvature remains well defined. The loop quantum gravity approach utilized here 
relies on densitized tetrads and SU(2) connections, and it is known that solutions with 
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classically degenerate metrics can yield finite equations when using these new Ashtekar 
variables [35]. Therefore, even before implementing quantum corrections it is worthwhile 
using these alternate variables as they are better suited to study systems with non-trivial 
topology. In fact, it is possible that degenerate tetrads may play an important role in 
quantum gravity [28], [29] . 




a) b) 

Figure 2: Schematics of topology changing space-times via wormhole formation. Figure a) 
represents topology change via the formation of an inter-universe wormhole. Fig- 
ure b) represents topology change via the formation of an intra-universe worm- 
hole. The points p represent the critical point of the topology change, where 
quantum gravitation effects are expected to become important. This sort of 
topology change is thought to be ubiquitous at the Planck scale. 



The cases for topology change or the existence of a space-time foam are still open. It 
is not yet completely clear if either of these are indeed possible or not. However, they 
are actually related issues as the existence of a space-time foam clS cl SG3i of fluctuating 
topologies is directly linked to the possibility of topology change. That is, if this picture 
of the space-time foam is correct, then topology change is allowed, at least in the quantum 
gravity realm. Although we also cannot answer this question definitively here, we do find 
that the existence of non-trivial topologies is less exotic when loop quantum corrections 
are present in the equations of motion, and hence this may be an indication that in the 
full loop-quantum gravity picture topology change would be at least more probable than 
in the classical regime. In Wheeler-DeWitt theory, a study has also been performed by 
Visser [30] to determine whether the existence of non-trivial topologies is favored or not in 
that theory. 

As well, the existence of wormholes in gravity theories which deviate from general rela- 
tivity in the high curvature (high energy) realm are of interest to the wormhole community. 
Of particular recent interest has been f{R) theory, where higher curvature scalar terms 
may appear in the gravitational Lagrangian. (See, for example, [31j.) Wormholes have 
been studied in the /(i?) theory [32j, [33j, [34J, [35j. The study in this article can be seen 
as another extension of general relativity to the high energy realm. Specifically related to 
discerning effects from quantum gravity, some studies of wormholes within the paradigm 
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of Wheeler-DeWitt theory have been carried out in [7], [l6]. As weh, many authors have 
used classical higher derivative theories of gravity to capture some possible effects of high 
energy corrections to the general relativity solutions [47J , [48j , [49j , some of which fall in the 
f{R) category mentioned above. Others still have approached the problem by considering 
semi-classical gravity where the non-trivial topologies are supported by quantum matter 

M - M- 

Given that the above solutions are necessarily anisotropic, we aim to study here in- 
homogeneous anisotropic systems. As most of the studies of such solutions have to date 
been classical, we also add some low-order loop quantum gravity holonomy corrections. 
We use a midi-superspace approach where we freeze the symmetry a priori to be spherical 
symmetry and then apply an "effective quantization" by replacing the configuration vari- 
able (the SU(2) connection) with a function incorporating some holonomy corrections from 
the quantum theory. Although we analyze anisotropic structures in general, we will pay 
special attention to the wormhole solutions, given the above arguments for their possible 
importance in quantum gravity. 

Regarding the possible detection of quantum gravity effects at the macroscopic level, 
it was suggested in Ref. [55] that one could imagine an extremely advanced civilization 
|56j pulling a wormhole from this submicroscopic space-time quantum foam and enlarging 
it to macroscopic dimensions. However, in a more plausible scenario, the possibility that 
inflation might provide a natural mechanism for the enlargement of such wormholes to 
macroscopic size was explored [57] . Some authors have also investigated the effects of such 
a foamy space on the cosmological constant. One example is the Coleman mechanism, 
where micro- wormhole contributions suppress the cosmological constant, explaining its 
small observed value [58]. However, even in the absence of detectable effects, predictions 
of a quantum gravity theory are important for understanding physics at the Planck scale 
[59]. 

II : SPHERICALLY SYMMETRIC ANISOTROPIC STRUCTURES AND 

QUANTUM CORRECTIONS 

The basic variables in loop quantum gravity (fixed to the time gauge) are the SU(2) con- 
nection as the configuration variable and the densitized triad as the conjugate momentum 
variable. In terms of more traditional variables in the 3-1-1 space-time decomposition, the 
densitized triad, E-"", is directly related to the spatial three- metric, Qab, via 

=det{q)q''^ . (1) 

In this work we adopt the convention that indices a, b, c, etc. denote spatial components 
whereas indices i, j, k, etc. are internal indices coupling to the su(2) algebra. 

The SU(2) connection encodes information about the extrinsic curvature of the spatial 
3-surface via 

A\:=r\ + jKi, (2) 

where 7 is the Barbero-Immirzi parameter. Here, T^^ is the spin connection whose associ- 
ated derivative annihilates the regular (i.e. non-densitized) inverse-triad, e*^: 

d^,e\^+e^^,T\^e%=0, (3) 
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and K'^a := . ^ 6^^ E J' Kab {Kab being the usual extrinsic curvature of the 3-surface). 
r*^ may be explicitly solved for via inversion of ^ as 



dae\-dbe\ + 5'''6^nei'eZdbe'', . (4) 



Writing the Einstein-Hilbert action in the 3+1 form of the ADM formalism in terms 
of A and E, and varying with respect to the shift- vector, A^**, and the lapse function, A^, 
yields the gravitational vector constraint and scalar constraint respectively: 



yj--) =E,''Fl, = 0, (5i) 
e'\Fi, - 2(1 + 7')K\kQ + 2AyditM| = 0, (5ii) 



^(grav) _ ) ' J 



yd^tOE) 

where we have included the cosmological constant (A) term for generality (which we set 
to zero in the subsequent analysis). The field-strength is given by F"^^^ := daA^^ — dbA^^ + 
e',f.A\A\. To this set one adds the Gauss constraint in order to eliminate redundant 
degrees of freedom from internal rotations: 

G'r^ := daE.^ + e,/A^,S,'^ =: D^Et = 0. (6) 

The formal quantization occurs by replacing the "classical functions" A*„ and E^ in 
([5il ii) and ([6]) with operators which obey the non-canonical algebra of loop quantum gravity. 
The resulting constraint operators then serve to constrain a ray in the Hilbert space with 
an appropriate inner-product [60], J6T], [62]. 

As dealing within the full quantum theory is technically difficult, we will deal with 
an effective quantum theory such as that which is often used in studies of loop quantum 
cosmology and the study of vacuum black hole interiors. That is, the connection shall be 
replaced with functions of A which can more easily be represented by holonomies. Details 
on this procedure shall be provided later. Specifically, it has been shown that, within 
this approximation, the true singularities present in the classical theory are alleviated in 
the corresponding quantum theory in both black hole interiors as well as in cosmological 
models. (See [63], |64] and references therein for details). As well, much progress has 
been made regarding the partial quantization of spherically symmetric inhomogeneous 
pure gravity in the interesting works on the exterior Schwarzschild space-time [65], [66] . 
[67], [68]. 

The systems studied here are necessarily inhomogeneous, anisotropic, and non-vacuum, 
which leads to complications not present in simpler systems. To begin with, we will require 
an appropriate inhomogeneous and anisotropic matter contribution to the equations ([5ilii) 
and ([6]) . For this we choose an anisotropic fluid whose four-dimensional stress-energy tensor 
is given by: 

T^u = (p + P±)Uf,Uu + p± g^v + {P\\ - P±)Sf_iSu . (7) 

Here p, p^ and p^ are the energy density, perpendicular (to the inhomogeneous direction) 
pressure and parallel pressure respectively as measured in the fluid element's rest frame. 
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The vector is the fluid 4- velocity and s'^ is a space-hke vector orthogonal to u^. On-shell 
these vectors satisfy: 

u% = -l, s% = +l, n% = 0. (8) 

The fluid is particularly useful in anisotropic studies as it provides one of the most 
general matter models allowed, yet still respecting the symmetry constraints. Essentially, 
for anisotropic structures in spherical symmetry we require a stress-energy tensor of Segre 
characteristic [1, 1, (1, 1)], and the spherically symmetric anisotropic fluid is one of the 
most general matter models which provide this. It may also be noted that the material ([7]) 
is also capable of accommodated Segre characteristic [1, (1, 1, 1)] (i.e. the perfect fluid) in 
the special limit p± = and can therefore be seen as a generalization on perfect fluids. 
A possibly more fundamental material, such as the scalar field, under certain assumptions 
on the field configuration which are compatible with spherical symmetry, is algebraically 
similar (in its stress-energy tensor) to the anisotropic or isotropic fluid [69j, and therefore 
by choosing a fluid source we potentially cover these scenarios as well. Finally, using the 
anisotropic fluid will simplify comparisons with purely classical solutions which most often 
use a fluid material as their source. 

Hamiltonians for a dust and an isotropic perfect fluid have already been constructed 
in the literature [70]- [76] and, as is usual in fluid mechanics, variables more appropriate 
for thermodynamic studies are usually chosen. However, here we choose to work in the 
original fluid variables of ^ as this system is more perspicuous and has immediate physical 
interpretation in these variables. The price to pay is that certain relationships between 
thermodynamic quantities (such as temperature, entropy, etc.) are not enforced from 
the variational principle in this scheme, but these are usually not of primary interest in 
gravitational studies. For an in-depth treatment of more fundamental held couplings to 
loop quantum gravity see [77], [75] . 

As a fluid of this type is considered an effective matter model, there seems to be no 
Lagrangian in the literature describing the anisotropic fluid. Therefore we adopt here 
a pragmatic strategy often used in held theories [79]. Specifically, we will construct a 
Lagrangian density, C^f, which produces the desired equations of motion, in this case the 
stress-energy tensor ([7]). After the derivation, the conservation law (which results from 
having a true scalar Lagrangian and the gravitational field equations) supplemented with 
the enforcement of the on-shell conditions ([8|) , which are also consequences of the variational 
principle in this scheme, will guarantee that the fiuid equations of motion are obeyed. We 
propose the following matter action which produces the desired result: 

■ [ I© 5"^ + ll - 2p^ + A - ll I d'^x , (9) 

Jm4 L J L J J 



=87r 

where we define G := (/? + pj_) and A := (py — p^), which must be viewed as independent 
variables in their own right (the quantity can be viewed as the particle number times the 
enthalpy). The quantities p and are true scalars as can be easily checked using ([7]) and 
([8]), via: T^^u^Uy = p and T^'^s^s^ = . Now, having established that p and are scalars, 
the trace T'^ = —p + 2p^ establishes that must also be a scalar quantity (and hence 
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so are and A). Note that with this particular choice of action the fluid variables and 
A appear simply as Lagrange multiplier fields which give rise to the set of constraints ([8]), 
and their dynamics are not completely governed by their variation. Instead, as mentioned 
above, the fiuid dynamics are governed by variation of the gravitational degrees of freedom 
via the Bianchi identities. This is in-line with the usual treatment of fluids in classical 
general relativity where one does not specify the fluid equations of motion separately. 

Next, all 4-dimensional quantities need to be split into a 3+1 decomposition, and the 
4-metric components replaced by functions of N, N"" and q°'^{Ef). For this we write 
[u^] = [u, Ua] and [s^] = [s, Sa]- After some calculation, the matter action in 3+1 form in 
the appropriate variables is found to be: 

UaUb - ^ + 1 



r (' r ATo / TP a TP b 



2px + A 



^2 J ^2 

^'''n^ + Vdi^M " ''''' " ^ " ' 



d^xdt. (10) 



Having now established a matter field in the correct variables, we freeze the symmetry 
of the system to spherical symmetry. For studies in the Ashtekar variables, we require an 
ansatz for a connection, A = A^^ Tidx"" , as well as a densitized triad, E = E^"" T^da , which 
is capable of accommodating this symmetry. We therefore utilize the following pair for 
this: 

A =AuiTi dy + {AiT2 + AiT-s) dO + (^^"3 - AuT2) sin(6') # + cos(6')ri d<j) , (Hi) 

d d d 

E =£uiri sin(0)— + {£,T2 + SuTs) sin(0)— + (^iTg - £:„r2) — , (llii) 
oy oO 0(p 

with 0<i;^<27r,0<^<7r and the functions A,, and are functions of the inhomoge- 
neous coordinate, y and time only. The Tj represent the standard su(2) generators. 

By utilizing (|llip and (jlliip in ([5i|) and ([6]) , and adding to them the matter contributions 
(from the variations of (10) with respect to the shift- vector for the vector constraint) we 
arrive at the following pair of constraints for our system: 

V =2sm{e) {Am [SiAu - £uA,] - £,A[ - £uA[,} 

+ levryl^ • |sin(6l)| ■ E ■ [Q U ■ + A ■ S ■ s^] = , (12i) 

G = sm{e) [2A,£u - 2Au£i + ^m] = , (12ii) 

where primes denote derivatives with respect to the inhomogeneous coordinate y, and we 
have dropped the index on V and G as each constraint only yields one equation. Here 
we have used the definitions E := ^J£^+£^^, U := u/N"^ and S := s/N'^ to simplify 
expressions. Also, we have set the shift vector to zero after variation as there is no loss 
of generality by doing this since we are restricting ourselves to spherical symmetry. Note 
that the fluid does not contribute to the Gauss constraint. Also of note is that in the case 
of a static fluid, the fluid makes no contribution to the vector constrain10. 



^It may be of interest to note that the staticity of the resulting space-time is actually encoded in this 
statement, although it is not obvious. Time derivatives of the metric are encoded in the vector constraint 
via the connection A^a through the extrinsic curvature in ((2]). If the matter contribution to the vector 
constraint vanishes, the only way to satisfy the spherically symmetric gravitational vector constraint is to 
make these terms involving the time derivative of the metric equal to zero. 
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For the scalar constraint the issue is more comphcated. This is because in the gravi- 
tational part of S, ()5iip . the extrinsic curvature terms appear and the extrinsic curvature 
needs to be replaced with the connection A via ([2]). Therefore we write 



(13) 



However, r*^ is a function of the triad and inverse triad via ([4]), and this must be ex- 
pressed in terms of the densitized triad E^^ only. Therefore we need to utilize the following 
relationships in (j4]): 



This yields the following spin-connection: 



TditM [E^-' ■ 



(14) 



1 F' f sinful 



2E^ 



+ cos(6')rid0 . 



(15) 



We now have all the quantities required to write the full scalar constraint in these 
variables as 

c 2|sin(6')| 



2£,£u,A[, - 2£u£uiA[ + 2£u,Au {£iA + £uAu) + (A^ - l) 

4e4a2 + e2 (£:;„)' + 4e2£:;„ {£,,a - £iAu) 

+ 8£uj (E^Aii + £u£i - S'Su) {£uAu + £jA,] 
+ Wtt^/\£^\ • |sin(0)| • E • [9 [/ • u + A • S • s - px] , 



(l + 7^)|sin(g)| 
2^/^-E372 



(16) 



where the definition A := y\Af+A[^ , and the fluid constraints dS]), have been used after 
variation. 

For the effective quantization we utilize a scheme that has proved fruitful in studies 
of loop quantum cosmology [63] and black holes [80], [81], [82]. That is, we replace the 
connection with functions that reflect the holonomy structure of the representation of the 
quantum algebra. The operator A is constructed from turning the Poisson brackets to 
commutators: 

A^,(x), ^/(y)] =ih6tS}6H^,y)- (17) 

Although j4*„ represented as "multiplication by can be used to satisfy the above 

bracket, it turns out that the operator A^^ does not have a well defined representation in 
the Hilbert space of connections. However, the space of functions of holonomies (cylindrical 
functions) does have a well defined, metric independent measure and can be turned into a 
mathematically rigorous Hilbert space. The holonomy of A is given by 



he{A) = V exp 



A 



(18) 
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where V represents a path ordering of the holonomy paths, e. In this vein we make the 
substitution, 

_^^^sm(A^^ (19) 
oj 

where 5j represents the length of the holonomy path. This holonomy path is not constant, 
but related to the proper area of the holonomy loop, as has been shown in studies of loop 
quantum cosmology and black holes [82], [83], [H]. More details of this will be provided 
later. Admittedly, this method is more rigorously motivated in certain studies of loop 
quantum cosmology [H], [85] . 

Although not required for the study of the constraints, for completion we construct the 
fluid Hamiltonian by identifying appropriate canonical momenta for the fluid: 



:= -167riVv^det(-E;) - Q-u^, := -167riVA/det(£;) • A • , (20) 

(1) (2) 

which leads, via the standard Legendre transformation, to the matter Hamiltonian 

i^af = j n^td'^x = 167r J Ny^\-\smi9)\ • E • [9 - A - p^] d^x . (21) 

(Care must be taken with minus signs as our Lagrangians are "minus" the more common 
convention.) In the case of the perfect fluid, the above reduces simply to an integral over 
the proper-energy density of the fluid. 

The above analysis is general, save for the assumption of spherical symmetry imposed 
from equations (jllip and (jlliip onward, and rather complex. We next simplify the system 
to study specific cases of interest. 



II. 1 : The static scenarios 

The models we shall study here will be time independent. Under the added assumption 
of staticity, the vector constraint reduces simply to the gravitational vector constraint, as 
a static fluid makes no contribution to the vector constraint. We shall also impose the 
following simplification on the connection and densitized triad, which is still compatible 
with spherical symmetry at least in the static case: 

A =AiiT3 dO - AuT2 sm{e) d4> + cos{9)ti dcj) , (22i) 

E =£uiTi sm{e)— + £,T2 sm{9)— + £,ts— , (22ii) 
ay oO o(p 

where now the functions A and £ are functions of y only. The relationship between the 
functions above and the usual metric in curvature coordinates is the following (with y in 
this case identified as the usual radial coordinate r ) from the relation det{q) ■ q""^ = Ef' ■ E^ : 

da^ = Qat dx^ dx^ =B{r) dr^ + dO^ + sin^ 6 dc/)"^ 

= dr"^ + £m de^ + £ni sin^ dcf)'^ . (23) 
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Prom this choice, the remaining gravitational vector constraint is identicahy satisfied and 
the Gauss constraint (jl2iip is satisfied by employing the following relation: 



(24) 



At this stage it is worth noting a few issues regarding the symplectic structure. By 
making the above simplification we are merely going to another set of canonical coordinates. 
To see this, one can make the following transformation on the original set of variables 



Ai =: cos(/3) , An =■■ A^ sin(/3) , 



(25) 



for some angle /3. Our particular choice corresponds to /3 = 7r/2. Similarly, via rotations 
on the densitized triad (both spatial and SU(2)), a canonical momentum can be defined 
via 



28^ =:2E'^ . 



(26) 



These configuration-momentum coordinates turn out to be exactly those of Bojowald and 
Swiderski [86], and Campiglia, Gambini and Pullin [66], [68] with symplectic structure: 



{^^(x),7r*(x')}=27<5(x 



X 



(27) 



Below we shall continue using the notation and for these components. 

Using the above simplification yields only a single connection component, ^n, that 
needs to be effectively quantized. Using the method discussed earlier, we make the substi- 
tution 

sin (Ai (^ii) 



5„ 



(28) 



where again (Jn represents the length of the holonomy path. As mentioned previously, it 
has been shown in a cosmological setting in [83] |84] and for black hole interiors in [82], 
that a reasonable semi-classical limit might not be obtained utilizing a fixed b in the above 
substitution. Instead, one can relate the length of the holonomy paths (the (5's) to the 
classical area. The proper area spanned by a holonomy "square-loop" in the — (j) sector 
is given by: 



(29) 



By setting the above area equal to the smallest area predicted by the spectrum of the area 
operator of loop quantum gravity, we arrive at what is known as the Ji' schem^ [83]. In 
full loop quantum gravity, the infinitesimal area operator is given by 



da 



1/2 



and has the following spectrum: 



da\^) = STi^eJjp{jp + l)\S 



(30) 



(31) 



^That is, the "mu-bar-prime" scheme. The symbol /i is often used instead of 5 for the length of the 
holonomy path in loop quantum cosmology. 
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with the Planck length. The integer p denotes which puncture (or area element) is under 
consideration, and jp can take on half-integer values which represent spins carried by the 
punctures. The covariant normal vectors to the surface element are denoted by and n^. 
Setting jp equal to one-half yields the smallest predicted area, which we denote as a^i^- 



Omin := \/3 47r7^; 



(32) 



so that, by (f29]) . we have 5u = y/a~/\/^^i ■ The results presented below are valid for 
any value of a^in, provided it is small. The Immirzi parameter, 7, is determined by some 
means such as black hole entropy calculations ^STj- ^96j. Putting everything together, the 
effective quantization amounts to the replacement 



sin ( Ai • V£^, 



(33) 



with Omin provided by ([32]) . 

At this stage, the only constraint which remains to be satisfied is the scalar constraint 
(I16p . With the imposition of the other constraints and staticity, the scalar constraint 
simplifies to 



+ levrv^^ • |sin(e)| • l^il • [e [/ • u - px] . 



Derivatives can be treated as usual continuous derivatives; specifically, for An: 



(34) 



sin I All 



n 

\ , / AiiJlT^ 
sm I , -I- cos ' 



A-,-, 



Aw <?jjj 
2 £111 



(35) 



Alternatively, one can initially treat the derivative as a finite-difference, as suggested in 
1681 for vacuum black holes: 



1 

- < 

e 



sm 



V^'iii 



X+e 



sm I 1 



\/'^min 



(36) 



where we use a constant discretization, e , for the derivative as the derivative is taken with 
respect to coordinate distance, not a holonomy path. 
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II. 1.1 : Wormhole 



The above analysis in usual curvature coordinates is suitable for many spherically- 
symmetric anisotropic structures. However, as mentioned in the introduction, of particular 
importance in studies of quantum gravity is the wormhole. In curvature coordinates a sub- 
manifold of the wormhole is generated by considering a curve with the correct properties, 
and then creating a surface of revolution from this profile curve. A sample profile curve, 
along with the corresponding surface of revolution is illustrated in figure El Note that in 
this chart two profile curves are required, one for the "top-part" of the wormhole (P+(r)) 
and one for the "bottom-part" of the wormhole (P_(r)). 




Figure 3: Wormhole profile curve, P±{r), in the 6 — tt/2 submanifold. The wormhole is 
generated via rotation about the x-axis (inset). 



In terms of the profile curve function, P±{r), the densitized triad components read: 

1 /2 

£, = r^/B^ = r{l + [drP±{r)f] , £,u = r\ (37) 

with r > ro (refer to figure). Immediately one can see a potential problem with the variables 
chosen when discussing the wormhole scenario. Note that at the wormhole throat (ro in 
figure [3]) the derivative of P{r) becomes infinite, and hence the densitized triad component 
£i is badly behaved at the throat. This malady affects the metric formulation as well. 
Hence, an alternative to the usual curvature coordinates will be used here. We use a 
system of coordinates introduced in [97], and utilized in [98] and [49], which is better 
behaved at the throat. Essentially, we rotate the entire chart by 7r/2 as shown in figure HI 
In this new coordinate system, the spherically symmetric anzats (jllih and (llliih still 
hold and the equations in this chart are only slightly more complicated than in the usual 
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Figure 4: Wormhole profile curve in the 9 = tt/2 submanifold using the rotated coordinate 
system. The profile function is given by r = Q{x) — P~^(x) and the radius of the 
throat is Qo- As before, the wormhole is generated via rotation about the a;-axis 
(inset). Only a single profile function Q{x) is needed now and the densitized 
triad components are finite at the throat radius. 



curvature coordinates. The line element is now 

da^ = dx^ =C{x) dx^ + Q^{x) dO^ + Q^{x) sin^ 6 dcj? 

(£ \^ 

=^^dx^ + £uide'^ + £uism^Od(t)^ . (38) 
iiii 

Here we have used the same simplification for A and E as we used previously. Note that 
the advantage of this chart is that the relation 

1 /2 

8, = Q{x)^/C{x) = Q{x) {l + [d^x)?} , £ui = [Q{x)f . (39) 

yields finite components at the throat as the derivative of Q{x) is zero, not infinite there. 
As well, the coordinate x may span -co < x < oo and hence it is now possible to cover the 
wormhole geometry with just a single well-behaved chart (if the poles of the two-spheres, 
6 = and 6 = n, are excluded). 

Again in the static case, the vector constraint is satisfied. The Gauss constraint is 
also again satisfied by imposing the relation (j24p . However, this condition is slightly more 
complicated in this chart due to the fact that the derivative of Sm contains the derivative 
of g(x). 

We do not present all the details here as the procedure is similar to the above. In short, 
the effective quantization is accomplished by the same substitution as previously done: 



A: ^ — ^ ^ , (40) 
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and the simpHfied scalar constraint reads as before ((31 

The most sahent features of a wormhole occur in the throat region. For example, 
in classical wormholes, the necessary violation of energy conditions in static wormholes 
occurs in the neighborhood of the throat. In fact, in the classical scenarios if the matter 
field falls off sufficiently fast, or is patched to a vacuum solution, the properties far from 
the throat do not differ greatly from similar systems with trivial topology. As well, it has 
been argued in ^99j that the global topology is too limited a tool to study wormholes, and 
a local geometric analysis near the throat is generally more useful in discerning interesting 
properties of wormholes. We therefore now concentrate on the near throat region and 
study the quantum gravity corrections in this vicinity. 

In the neighborhood of the throat (x = 0), the profile function, Q{x) must have the 
following properties: 

i) Qo ■■= Q{0) > 0, 

ii) (5'(x)|^.=o = 0, 

iii) Q"{x) > in some neighborhood of the throalH. 

Aside from the above properties we make the mild assumption that Q{x) is analytic. Now, 
8i and £"111 are functions of Q{x) and its derivative via ([39|) and one can substitute these 
functions into (134^ to study the remaining scalar constraint. We shall first concentrate 
on the purely classical scenario, and hence simply replace in (j3^ by using (HH); i.e. 
= '^•iii/(2<?i)- This yields the following scalar constraint: 



_ 2|sin(g)| 



1 + (QO 



2Q-Q"-{Q'f-l +167r|sin(e)|-gyi + (Q0V- (41) 



We treat p as the unknown and hence the constraint 5 = yields: 



87rp 



[1 - 2QQ" + (Q- 

g2 [1 + {Q'\ 



,212 



(42) 



The above expression agrees with the expression one may derive by calculating the t — t 
component of the Einstein tensor using metric (j38p and ()39p using an arbitrary lapse (which 
does not appear in Next, the analyticity of Q{x) allows for a Taylor expansion about 
a; = 0, resulting in the following near-throat expression: 



1 



2Qi 



Qo 



1 r 



Qo + Qo"Qo - 4 (Q'oT Q, 



+ 



3Q^ 



Qo'QoQo - Qo" - QrQo + 24Q[;' (Q'o')' Qo + 0{x 



(43) 



^More precisely, if Q's first non-zero derivative (higher than first order) at x = is of even order, the 
function attains a local minimum if this derivative is positive, and hence we have a wormhole throat. If 
its first non-zero derivative is of odd order, it is a point of infiection and therefore does not describe a 
wormhole throat. 
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where the subscript indicates that the quantity is evaluated at x = 0. (We do not demand 
that terms odd in powers of x vanish as the wormhole does not need to be symmetric about 
the throat.) 

Next we tackle the more difficult case with the quantum holonomy corrections. In this 
case, in is first replaced via (|10|) . and the derivative of is replaced by (j35]) . 



after which the substitution An = £l^^/(2£i) (from the Gauss constraint, ([MD) is utilized. 
The full analytic expression for p can be calculated but the result is quite complicated so 
we simply present the series result: 



Sfrp 



\ (1 - 2Q'o'Qo) - 



+ 







1 



QO + Qo'Ql - 4 (Qo)' - «n.in (Qo 



X 



1 



3Q8 



(44) 



The purely classical limit is achieved when a^^^ — ?• 0, which agrees with ()43p . Note that the 
quantum corrections (terms multiplied by 0^1,1) do not contribute exactly at the wormhole 
throat. Instead, they come in at order and, given the sign of the second derivative of 
Q{x) near the throat, the quantum corrections contribute positively to the energy density 
near the throat and hence may act to lessen energy condition violation for the matter in the 
vicinity of the throat. Even though the energy density can be positive at the throat even 
in the purely classical scenario, it is known that somewhere in the vicinity of the throat 
energy conditions must be violated (see, for example, [99]). The fact that this modified 
theory of gravity increases the energy density is a positive indication for the lessening of 
energy condition violation, although the pressures should be also be studied for a concrete 
statement. We make some comments on the pressures below. 

Another expansion, although perhaps not as useful, is an expansion of p in powers of 
Oniin. This yields an expression that is order- by-order in powers of quantum corrections: 



Svrp 



[1 - 2QQ" + iQ') 



+ 



Q2[1 + (Q')2]2 3[l + (g/)2]3g4 



{Q'f + {Q'f - ^Q"Q 



I8O72 [l + (Q')2]4g6L 

{Q'f 

252O72 [1 + (Q')2]5q8 



(97^ + 5) {{Q'f + {Q'f) - i5i'Q"Q 



^min 



mm 



(57' + 7) {{Q'f + {Q'f) - 7j'Q"q] aL + 0{atJ 



(45) 



Compatible with the previous result, it can be noted that at local extrema where Q' = 
(eg. the throat), terms arising from the quantum holonomy corrections do not make a 
contribution, but do contribute slightly away from the extremal point. It is therefore 
possible in principle to have a model where the energy densities of both the classical and 
quantum corrected throat are zero at the throat, but in the vicinity of the throat region 
the energy density of the classical solution is negative and the quantum one is positive. 
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Ideally one would like to analyze the individual pressures to completely specify the 
properties of the matter field. However, in the 3+1 connection formalism this is not readily 
available^. In the case of isotropic loop quantum cosmology with a scalar field, one can 
extract information on the (isotropic) pressure and this has been achieved in |100j where 
energy conditions were studied. 

Finally we present a few specific models to compare the quantum corrected models 
with the purely classical models. In figures [5^-d) a specific wormhole profile Q{x) = 
ao cosh(x/2;o) is chosen and the quantum corrected (solid lines) and classical (dashed) 
energy densities are plotted for various values of the parameters, with both positive and 
negative energy densities near the throat. Note that the quantum effects tend to raise 
the energy of the fluid, and hence the exotic nature of the fluid tends to be lessened in 
comparison to the purely classical scenarios in the cases where the energy densities are 
negative. In figures [6^-b) we also plot non-symmetric wormholes with a profile given by 
Q{x) = ao cosh(x/a;o)+/3o2;^ for sufficiently small (3q so as not to spoil the local minimum. It 
can also be seen here that the quantum corrected versions tend to have more positive energy 
densities. In all cases where negative energy densities occur, the region of negative energy 
is smaller in the corresponding quantum corrected case. It is possible that those particular 
models which have outwardly increasing energy density may be unstable. However, if these 
models are to represent ephemeral quantum fiuctuations in the vacuum, it may be that 
instability is not a serious issue. 
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Figure 5: Symmetric wormhole models with Q{x) = ao cosh(x/a;o). The parameters are as 
follows: Figure a): ao = 0.38, xq = 0.86, figure b): A close up of the throat region 
of the previous figure, figure c): ag = 0.005, xq = 0.06, d): ao — 0.3537, xq = 0.5. 
In all cases 7 « 0.27 and ip was set to 0.1 to exaggerate the differences to make 
them easier to see. 



In classical general relativity, there are singularity theorems which state that if the clas- 
sical equations of motion hold in scenarios describing objects such as gravitational collapse 

*This would involve writing the 3+1 action (including the connection terms) completely in terms of the 
metric and its derivatives, vary the action with respect to the metric, then re-write terms which can be 
expressed as connections again as connection terms, so that the substitution (|33|) may be performed to get 
the quantum corrected pressures. 
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I ^ \l \ 

a) b) 

Figure 6: Non-symmetric wormhole models with Q{x) = ao cosh(a;/a;o) + Pqx^ . The pa- 
rameters are as follows: Figure a): ao — 0.25, xq = 0.5, /3o = 0.5, figure b): 
ao = 0.305, xo = 0.4, /3o = —0.95. In both cases 7 ~ 0.27 and was set to O.f 
to exaggerate the differences to make them easier to see. 



to black holes or collapsing cosmologies, and the material present obeys energy conditions, 
then a singularity is inevitable [lOlj . |102j . One loop-hole out of these arguments is the 
abandoning of the energy conditions, as it is known that energy condition violating mat- 
ter can theoretically prevent such singularities from forming. It is also often believed that 
quantum gravity effects should take over close to singularity formation, when the curvature 
is large, and that these effects may prevent the formation of the classical singularity. In this 
vein, the loop quantum gravity paradigm has previously indicated that both in cosmology 
and black holes the quantum gravity effects have replaced the singularity with a smooth 
bounce [T03j . [104j . [80j . [82j due to gravity becoming repulsive under certain situations in 
the quantum regime. That is, a scenario that classically requires matter energy condition 
violation to occur, can occur naturally within loop quantum gravity. What we find here in 
the context of wormholes is a manifestation of the same phenomenon. In order to support 
the wormhole classically, energy conditions must be violated to a certain degree. In the 
quantum corrected case, however, the energy condition violation is less, and this is directly 
due to the quantum properties of the gravitational field. It is, in principle from the above 
expressions, possible to find wormhole solutions where there is an energy density which 
is everywhere positive (with a zero at the throat) but whose classical counter-part yields 
an everywhere negative (with a zero at the throat) fluid energy density. This is due the 
repulsive nature of gravity in the regime where loop quantum effects become important. It 
should be noted that these are low-order quantum corrections, similar to those employed 
in most studies of effective loop quantum black holes or loop quantum cosmology. It is a 
possibility that higher order corrections may improve the energy condition situation even 
further. Related to this, the repulsive nature of the gravitation under such extreme con- 
ditions could act to alleviate the singularities that are thought to occur when space-time 
topology changes classically. (Although even classically it has been shown that there are 
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some topology changing scenarios where the singularities are not curvature singularities 
but merely metric ones j27j.) 

Ill : CONCLUDING REMARKS 

We have studied anisotropic spherically symmetric systems in the 3+1 Hamiltonian for- 
malism of gravity in the variables used in loop quantum gravity. By replacing the SU(2) 
connection with functions that encode the holonomy structure of the corresponding oper- 
ator we essentially have an effective theory with some quantum inspired corrections. In 
the case of wormhole throats it was found that the energy density of the material source 
is increased in comparison to the purely classical case. This may indicate that the energy 
condition violation ubiquitous in wormhole throats in Einstein gravity may be lessened 
by quantum gravity effects. This provides another arena to study higher energy effects 
of gravity, complementing other techniques such as those incorporating higher curvature 
effects through a modified gravitational Lagrangian (such as f{R) theories [39], |31j-|35j). 
As well, the wormhole provides yet another theater to test the predictions of a quantum 
gravity theory. 
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